Multiple-gap structure in electric-field-induced surface superconductivity 
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Local superconducting gap structure is studied as a function of nanoscale depth in electric-field- 
induced surface superconductivity such as in SrTiOa. We examine solutions of Bogoliubov-de Gennes 
equation in two limiting confinement potential cases of electric field with and without screening 
effects. As unique properties different from bulk superconductivity, there appear in-gap states 
even for isotropic s-wave pairing, due to multiple gap structure of sub-band dependent surface 
superconductivity. These determine the depth-dependence of local superconductivity. 
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I. INTRODUCTION 

Electric-field-induced carrier-doping technique, using 
field-effect-transistor (FET) structure or electric-doublc- 
layer-transistor (EDLT) structurefS^ attracts much at- 
tention as a new method to carrier doping, other than 
the methods of chemical doping. A merit of electric-field- 
induced doping is that we can control the doping carrier 
density by a gate voltage in a same sample. This will be a 
powerful platform in future studies of condensed matter 
physics. When this is used at the surface of insulators, 
carriers are induced near the surface, and trapped in the 
confinement potential of the electric field. Using enough 
strong field by EDLT, we can realize superconductiv- 
ity of the surface metallic states at low temperature T, 
such as in SrTiOa, 111 ZrNClplKTaOaplMoSa!^ The gate 
voltage control of surface superconductivity in SrTiOa 
was also realized at the interface of LaTi03 /SrTiOa and 
LaAlOa/SrTiOa. 1113 

Compared to these developments of experimental re- 
search, theoretical understandings are not enough for 
properties of the electric-field-induced surface supercon- 
ductivity. We have to discuss whether the surface super- 
conductivity has the same properties to those of bulk su- 
perconductivity, or whether it has different unique prop- 
erties. In future experiments for unconventional super- 
conductivity produced on a surface, we have to distin- 
guish unique properties of surface superconductivity and 
exotic properties of unconventional superconductivity. 

As different properties from bulk metallic states, sub- 
bands are formed in the surface metallic states due to 
the confinement potential by strong electric fields.^ Since 
multiple sub-bands are occupied by surface carriers, this 
system is not ideal two-dimensional states. We also note 
that the local carrier density n(z) has spatial variation 
as a function of depth z from the surface in the surface 
metallic states, while n(z) is constant in bulk metallic 
states. The quantitative estimate for the z dependence 
is one of problems for the electric-field-induced metallic 
state. Therefore, also in the theoretical studies of sur- 
face superconductivity, we need to know detailed spatial 
structure of the superconducting gap in the nanoscale 
and its sub-band dependence. These studies enable us to 
find differences from bulk superconductivity. 



In this letter, we study unique properties; local elec- 
tronic states and sub-band dependence in electric-field- 
induced superconductors. We will discuss multiple-gap 
structure of the sub-band dependent surface supercon- 
ductivity. Since we determine the spatial structure in the 
order of Thomas-Fermi length near the surface, we solve 
the Bogoliubov-de Gennes (BdG) equation^ under the 
electric-field F(z). We discuss the depth z dependence 
perpendicular to the surface at z = 0. As for confinement 
potential V(z) by F(z), we compare two case s; triangu- 
lar potential and self-consistent potentialP^ The latter 
is the case when induced carriers completely screen the 
applied electric field. The former is the opposite limit 
where the screening is negligible. 

This paper is organized as follows. After we explain 
our theoretical formulation of BdG equation under elec- 
tric fields in Sec. [TlJ we study the depth-dependence 
of local superconducting gap structure in Sec. |III| and 
the gap structure of the sub-band modes in the spec- 
tral weight in Sec. 



IV 



In order to discuss the rela- 
tion of sub-band dependent gap structure and the depth- 
dependence of superconducting states, we perform the 
analyses of sub-band decomposition for the surface su- 
perconductivity in Sec. |Vj The last section devotes 
to discussion and summary, including the topics of the 
Bardeen-Cooper-Schrieffer(BCS) - Bose Einstein conden- 
sation(BEC) crossover phenomena in the surface super- 
conductivity. 



II. BOGOLIUBOV-DE GENNES THEORY IN 
CONFINEMENT BY ELECTRIC FIELD 

Throughout this letter, energy, length, local carrier 
densities are, respectively, presented in unit of eV, nm, 
and nm -3 . We typically consider the case of sheet carrier 
density n 2 D = 6.5 x 10 13 [cm -2 ], and electric field at the 
surface is given by F Q = F(z = 0) = 1.4 x 10~ 3 [V/nm]. 
The triangular potential with this Fq corresponds to one 
of the case calculated in Ref. [1] for SrTiOa , while single 
band case of effective mass to* = 4.8too is considered 
here, Too is f ree electron's mass. 

In the normal statep^ the eigen-energy E t and wave 
function u c (r) = e^ kaiX+kvV ^u e (z)/^/S are determined by 



the Schrodinger equation 

Ku e (z) = E e u e {z), 

with 



h 2 d 2 

K = —_ ^+Eu+V(z)-n, 



2m* dz 2 



(1) 



(2) 



k 2 + k 2 , and S is unit area of 



E\\ = Ti 2 k 2 /2m*, fcg - i- 2 ->- ^ 
surface. We assume that the wave functions vanish at 
z = as the boundary condition. In the parallel direction 
to the surface, the eigen states are given by wave numbers 
k x and k y of plane waves. Thus, the eigen-states of Eq. 
(Til) are labeled by e = (k x ,k y ,e z ). e z indicates label for 
sub-bands coming from quantization by confinement in 
the z-direction. The local carrier density is calculated as 
n(z) = 2n T {z) with 



(z)=5> £ (z)| 2 /(£ £ ), 



(3) 



where f(E) is the Fermi distribution function. To fix 
n 2Bi we tune chemical potential /j,. In the triangu- 
lar potential case, the confinement potential is given by 
V(z) = \e\FoZ. For the selfconsistent potential, 

F(z) = F (l - J n(z')dz'/n 2 n) , (4) 

by Gauss's law, considering the screening by n(z), and 

V(z) = |e| [ F(z')dz'. (5) 



As F(z — > oo) = 0, ri2D = Jo n(z')dz' . Iterating cal- 
culations of Eqs. ([TJ)- ^ and Eqs.(|3])-([5| in the region 
< z < L, we determine V(z) in the case of selfcon- 
sistent potential. Typically we use L = 80[nm]. 
In the superconducting state, the wave function 



u e (r) 
i> e (r) 



1 



i(k m x+k y y) ( U e (z) 



is determined by solving the BdG equatio n 1 8 ! 11 ! 12 ! 



K A(z) 
A(z) -K 



u e (z) 
v e (z) 



= E f 



u t (z) 
v e (z) 



(6) 



(7) 



The pair potential A(z) is selfconsistently calculated by 
A(z) = ^ pair J2 u € (z)v e (z)f(-E e ) (8) 



with the energy cutoff E cut of the pairing interaction. 
Here, we consider a conventional case of isotropic s-wave 
pairing. We typically use V^> a i r = 0.04, E cut — 0.01, and 
T~0. 
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FIG. 1: (Color online) (a) Density plot of LDOS N(E, z) as 
a function of z and E + fj,. Solid line presents the confinement 
potential V(z). (b) LDOS N(E, z) in (a) is focused near the 
superconducting gap. (c) DOS N(E) as a function of E + fj,. 
Dashed lines present N(E) for normal state. Left (Right) 
panels are for triangular potential (selfconsistent potential). 



III. DEPTH-DEPENDENCE OF LOCAL 
SUPERCONDUCTING GAP STRUCTURE 

First, we study the l ocal d ensity of states (LDOS) 
N(E, z) = 2N t (E, z) witfPE2l 



N t (E,z)= J2 



\u e (z)\ 2 S(E-E e ). (9) 



The left panel of Fig. [f|a) presents N(E, z) for the 
triangular potential. There, we see steps of LDOS by 
the sub-band structure of quantized bound states, as 
in the normal stated The lowest sub-band appears at 
E > Ei—x tW ia ~ 0.0059 near the surface. The continuum 
distribution above -Efci.min comes from finite E» > 0. 
Similarly there appears the LDOS of second sub-band at 
E > J5j= 2 ,min ~ 0.0103, and the LDOS of third sub-band 
at E > i?; = 3,min ~ 0.0138. Their contributions are over- 
lapped each other at higher energies. When the sub-band 
level I is higher, the eigen-energy Ei min becomes higher, 
and the distribution spread until deeper z from the sur- 
face. The superconducting gap appears near fj, ~ 0.0198. 

The N(E, z) for the self-consistent potential is pre- 
sented in the right-panel of Fig. [IJa). There, we see 
step-structures of sub-bands at low energy region, but 
the step size becomes smaller at higher energy, because 
the slope of V(z) decreases to zero as a function of z by 
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ity. For the sub-band decomposition, we calculate the 
spectral weight N(E, fcii , z) — 2Nf(E, fey , z) given by 



FIG. 2: (Color online) Density plot of spectral weight 
N(E, fey) as a function of E\\ = h 2 k^/2m* and E + fi for (a) 
triangular potential and (b) selfconsistent potential. Numbers 
in the figure indicate sub-band level I. 



the screening effect. Since the chemical potential is lo- 
cated at /i ~ V(z —> oo) ~ 0.0154, occupied states with 
E < fi are bound states, and empty states with E > fi 
are scattering states which are free from the confinement 
potential. The superconducting gap opens between the 
bound states and the scattering states. 

The superconducting gap structures are focused in Fig. 
|TJb) . Even in the isotropic s-wave pairing, we see in-gap 
states which have oscillations as a function of z and steps 
of gap-edges as a function of E, as characteristic features 
of electric-field-induced surface superconductivity. High 
intensity peaks of N(E, z) correspond to the maximum 
gap-edge, whose gap amplitude decreases discontinuously 
with increasing z. In the selfconsistent potential (right- 
panel), it reduces to zero at large z. 

Figure [TJc) shows density of states (DOS) N(E) after 
z-integration of N(E,z). Because of the in-gap states, 
gap structure in N(E) is different from that of bulk 
isotropic s-wave superconductors. In the triangular po- 
tential (left panel), the gap-edge has width from min- 
imum gap to gap-edge peak of maximum gap, as in 
anisotropic s-wave superconductors. In the selfconsistent 
potential (right panel), full-gap structure does not exist, 
since low energy states exist until near [i. The gap shape 
is similar to that of anisotropic superconductors with gap 
nodes. 



IV. GAP STRUCTURE OF SUB-BAND MODE 
IN SPECTRAL WEIGHT 

We discuss that these superconducting gap structures 
come from the sub-band dependence of superconductiv- 



N t (E,k l{ ,z) = ^Mz)| 2 <5(£-£ £ 



(10) 



from N t (E,k h z) 
function 



-it 1 lmG^(E,k\\,z) with Green's 



G t (E,ku,z)= / G ^ k ^ +k ^G^E,v,v')Axdy\ z=z (l\) 



G t (E,r,r')=J2 



E + iO - E e 



(12) 



and (x,y) = [x — x,y — y))^ The z-integration of 
N(E, fc|| , z) is given by N{E,k\\) = N(E, fey , z)dz. In 
Fig. [2]we show N(E, fey), which appear at the eigen en- 
ergies E e . There we see multiple parallel lines of the dis- 
persion relation as a function of E\\ , corresponding sub- 
bands of surface bound states. From the bottom, the 
lines are assigned to sub-band level I = 1, 2, • • •, as indi- 
cated in Fig. [2] In the case of self-consistent potential, 
the energy distance of the dispersion relation between 
sub-bands decreases for higher sub-bands, and the spec- 
tral weight becomes continuous near En ~ at E > fi in 
the scattering state. 

In the superconducting state, gaps open at crossing 
points of the particle mode and the inverted hole mode 
at E = /i, forming Bogoliubov's dispersion relations of 
superconductivity for each level of sub-band. The super- 
conducting gap is larger for lower sub-bands, indicating 
multiple-gap structure of surface superconductivity. In 
the case of triangular potential, the occupied lower sub- 
bands have different but finite gaps. In the case of self- 
consistent potential, sub-bands are occupied until quite 
higher levels, where superconducting gap reduces to zero. 
Therefore, the in-gap states appear until near E = /i. 
As for the z-dependence of spectral weight N(E, fey , z), 
the contribution of lower sub-band is dominant near the 
surface. The contributions of higher sub-bands become 
dominant at deeper z. 



V. SUB-BAND DECOMPOSITION OF LOCAL 
SUPERCONDUCTING STATES 

In Fig. [3]ja), we present the local carrier density n{z) 
and the sub-band decomposition. The eigen-states on 
the dispersion relations in Fig. [2] are classified to each 
sub-band level I. In Z-th sub-band, the wave function 
of the form of Airy functions has I — 1 nodes along z 
direction.^ The higher sub-band contributions can pene- 
trate into deeper z. Since the LDOS is integrated over 
-E/,min < E < jj, to obtain n(z), lower sub-band contribu- 
tions to n(z) becomes larger, because of smaller Ei tm - ln . 
The pair potential A(z) and the sub-band decomposi- 
tion in Fig. [3^b) have similar spatial structure to those 
of n(z). It is noted that sub-band dependent pair po- 
tential becomes smaller for higher sub-bands. In the 
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FIG. 3: (Color online) (a) Local carrier density n(z) and (b) 
superconducting pair potential A(z) as a function of z. The 
sub-band decompositions are also presented for sub-band lev- 
els I — 1, ■ • ■ , 5. Left (Right) panels are for triangular poten- 
tial (selfconsistent potential). 



sclfconsistent potential (right panels), while lower sub- 
band contributions are dominant, n(z) and A(z) includes 
contributions from further higher sub-band levels I > 5. 
Therefore, tails of n(z) and A(z) survive until deeper z. 

To discuss the origin of superconducting gap structure 
in Fig. [TJ sub-band decompositions of the LDOS are pre- 
sented in Fig. [4] From Fig. Qa), we see that lowest sub- 
band contribution (I = 1) to N(E, z) has large constant 
superconducting gap, but its distribution is restricted in 
very near the surface. The contributions from higher 
sub-bands I = 2, 3, • • • have smaller constant gap, and 
the distributions spread until deeper z. By the combi- 
nation of these sub-band contributions, the in-gap states 
and the z-dependence of superconducting gap structure 
in Fig. [l] are created. These sub-band contributions are 
clear also in the LDOS spectra in Figs. [4jb) and[4|c). 
There, we see multiple peaks of gap edge from sub-band 
contributions. Near the surface (z — 1.6), all sub-band 
contributions appear, and lower sub-band contributions 
are dominant. Therefore the main peak corresponds to 
gap edge of largest gap by the lowest sub-band. In lower 
panels of Figs. |4|b) and[4|c) for deeper z, since lower sub- 
band contributions (/ = 1 and 2) vanish, the main peak 
of gap edge appears at lower gap energy corresponding 
to higher sub-band (I = 3) contributions. 

In addition to the superconducting gap at E = \i, there 
appear extra small gaps at energies outside of supercon- 
ducting gap, as seen in Fig. [2j This occurs by the cross- 
ing of the hole- and particle-modes between different sub- 
bands. We see these extra gaps in higher sub-band con- 
tributions also in Fig. [41(a) . Because of the extra gaps, 
the LDOS in Figs, gb) and gc) has many extra peaks 
outside of the main gap energy. 

We note that low-energy in-gap electronic states are 
not determined locally by A(z) in the length order of 
nano-scale in this system. This is contrasted to conven- 
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FIG. 4: (Color online) (a) Density plot of sub-band decom- 
position of LDOS as a function of z and E + [i for sub-band 
levels 1=1, 2, ■ ■ •, 5 in the case of self-consistent potential, (b) 
LDOS N(E,z) and the sub-band decomposition at z = 1.6 
and 9 in the case of triangular potential, (c) The same as (b), 
but at z = 1.6 and 12 in the case of selfconsistent potential. 



tional case when A(z) is suppressed in the length scale of 
superconducting coherence lengthP^E^I There low-energy 
states appear as localized states by the suppression of 
the local gap. In the system of the electric field-induced 
surface superconductivity, approaching z — > near the 
surface, A(z) is suppressed toward zero in the length or- 
der of nano-scale, as shown in Fig. [3ja). However, from 
Fig. |4|a), we see that localized low energy in-gap states 
do not appear at the surface region (z < 2) of suppressed 
A(z). The local state of lowest sub-band has uniform 
gap with largest gap amplitude [top panel in Fig. [4|a)]. 
This indicates that the in-gap states reported in this pa- 
per is not due to the suppression A(z — > 0) — > 0. Rather, 
the in-gap states comes from the deeper region, as tails 
of wave functions for higher sub-band levels in Fig. 4][a). 
This is one of intrinsic natures in the electric field-induced 
superconductivity. 
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VI. DISCUSSION AND SUMMARY 

As future experiments to confirm the in-gap states due 
to the characteristic multiple gap structure, we expect 
observations of LDOS such as by point contact tunneling 
spectroscopy, which will see the gap structure in upper 
panels of Figs. 4(b) and 4(c). The contributions of the 
in-gap states will be observed in experiments sensitive 
to DOS of the superconducting gap structure, such as 
magnetic resonance, optical absorptions, etc. Electric- 
field-induced doping will be important platform to study 
unconventional superconductivity. Before that, it is im- 
portant to clarify the difference of properties between 
surface superconductivity and bulk superconductivity in 
conventional superconductors, as suggested in this work. 
As a concept of multiple-gap structure, electric-field- 
induced surface superconductivity can be said a new type 
of multi-band superconductors. We will see some similar 
behaviors to those of multi-band superconductors such 
as in MgB 2 and Fe-based superconductors. The number 
of contributing sub-bands can be controlled by the gate 
voltage. 

We point out an interesting pos sibility to realize the 
BCS-BEC crossover phenomena^fl^oi by controlling the 
gate voltage in surface superconductivity. In the cold 
atomic gases, the BCS-BEC crossover is see n by tun- 
ing the interaction via a Feshbach resonance! 21 * 22 ! The 
BCS-BEC crossover in multi-band superconductor was 
suggested by ARPES experiment in FeSe x Tei_ x PS The 



same situation appears in the surface superconductivity. 
In Fig. [2ja), the superconducting gap in 5th sub-band 
opens at the bottom of the band dispersion. That is, 
since the gap amplitude |A| is larger than the Fermi en- 
ergy Ep(= A* — -Ei=5,min) from the band bottom, the BEC 
regime |A| > Ep is realized. The gaps in other 1-4 th 
bands are in the BCS regime |A| < Ep. As mentioned 
above, the gap magnitude can be tuned by the gate volt- 
age. 

In summary, local superconducting gap structure and 
the sub-band dependence in electric-field-induced surface 
superconductivity were studied by solving microscopic 
BdG equation. There, the in-gap states appears due to 
the multiple gap structure of multiple sub-band super- 
conductivity, even for isotropic s-wave pairing. We evalu- 
ated how these structures depend on the screening condi- 
tion, i.e., triangular potential or selfconsistent potential. 
These different characters from the bulk superconduc- 
tivity, by the sub-band dependent multi-gap nature, are 
important to be considered, when we discuss properties 
of electric-ficld-induced surface superconductivity. 
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